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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  describes 
procedures  to  calculate  temperature  change  within  the  Coastal  Modeling  System  (CMS)  operated 
in  the  Surface-water  Modeling  System  (SMS),  version  12.1  (Aquaveo  2010).  The  defined 
procedures  are  demonstrated  in  an  application  of  temperature  modeling  in  the  Corrotoman  River, 
Virginia. 

INTRODUCTION:  Water  temperature  is  an  important  physical  property  influencing 
environmental  conditions  in  aquatic  systems.  Its  value  changes  with  atmospheric  conditions, 
typically  ranging  from  close  to  0  °C  in  the  polar  regions  to  more  than  30  °C  in  tropical  regions.  In 
coastal  zones  and  estuaries,  both  temporal  and  spatial  variations  in  water  temperature  are  also 
controlled  by  changes  in  circulation,  tide,  wind,  and  freshwater  inflows.  These  changes  in 
temperature  can  have  major  effects  on  water  density  and  water  stratification,  which  can  change 
ocean  mixing  and  modify  circulation  patterns.  Dynamic  behavior  of  suspended  sediment  can  be 
controlled  by  the  density-driven  flow  and  mixing.  Temperature  can  alter  the  water  physical 
environment  that  impacts  marine  organisms  with  the  change  of  water  turbidity  in  coastal  and 
estuarine  systems.  Water  quality  and  ecological  models  often  require  input  of  water  temperature 
and  salinity  information  from  a  hydrodynamic  model. 

COASTAL  MODELING  SYSTEM  (CMS):  The  CMS  calculates  water  levels,  currents,  and 
waves  through  the  coupling  between  a  hydrodynamic  model  (CMS-Flow)  and  a  wave  model 
(CMS-Wave).  These  two  models  can  be  coupled  to  simulate  sediment  transport,  morphology 
change,  and  salinity  transport  (Lin  et  al.  2008;  Sanchez  et  al.  2014;  Li  et  al.  2012). 

CMS-Flow  is  a  two-dimensional  (2D),  finite-volume  model  that  solves  the  mass  conservation 
and  momentum  equations  of  shallow  water  motion.  CMS-Flow  is  forced  by  water  surface 
elevation  (WSE)  (e.g.,  from  tide)  and  river  discharge  at  model  boundaries,  wave  radiation  stress, 
and  wind  forcing  over  a  model  computational  domain.  Physical  processes  calculated  by 
CMS-Flow  include  wave-current  interaction,  sediment  transport,  morphology  change,  and 
salinity  transport.  CMS-Wave  is  a  2D,  quasi-steady  spectral  wave  transformation  model.  The 
model  contains  theoretically  derived  approximations  of  wave  diffraction,  reflection,  and  wave- 
current  interactions  for  wave  simulations  at  coastal  inlets  with  jetties  and  breakwaters. 

TEMPERATURE  CALCULATIONS  IN  CMS:  The  CMS  calculates  the  temperature  field 
based  on  the  following  2D  heat  transfer  equation: 
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where  T  is  the  depth-averaged  temperature;  d  is  the  total  water  depth;  qx  and  qy  are  the  flow  per 
unit  width  in  the  x  and  y  directions,  respectively;  Kx  and  Kv  are  the  diffusion  or  mixing 
coefficients  of  heat  in  the  x  and  y  directions,  respectively;  p  is  the  water  density;  cp  is  the 
specific  heat;  and  Jt  is  the  net  heat  flux  across  water  or  bed  surfaces  into  the  water  body. 

Equation  (1)  represents  the  horizontal  heat  transfer  (advection  and  diffusion)  in  a  water  body 
considering  the  heat  exchanges  across  the  water  and  bed  surfaces  as  heat  sources/sinks.  The  heat 
content  in  the  water  body  is  primarily  controlled  by  air-water  heat  exchanges  whereas  the  heat 
flux  due  to  conduction  and  convection  across  the  bed  surface  is  assumed  to  be  negligible  and  is 
not  included  in  this  calculation.  In  addition,  heat  inflow  from  rivers  and  by  tidal  currents  also 
affect  the  temperature  in  the  water  body,  which  can  be  specified  as  boundary  conditions.  Heat 
flux  components  across  the  water  surface,  including  short-wave  solar  radiation,  long-wave 
atmospheric  radiation,  latent  heat  flux,  and  sensible  heat  flux,  can  be  calculated  in  the  CMS 
using  the  following  meteorological  parameters:  solar  radiation,  cloud  cover,  air  temperature, 
wind  speed,  and  surface  water  temperature.  The  bulk  formulas  to  evaluate  the  heat  flux 
components  were  described  by  Wu  (2008),  as  outlined  in  the  following: 

Short  Wave  Solar  Radiation,  Jsw 

J SW  =  J SW,CLR  (1  "  0.65C CLD X^-  _  ^ SW  X^  —  f SHD  )  (2) 

where  JswaR  is  the  short-wave  solar  radiation  received  on  a  water  surface  during  a  clear  (i.e., 
cloudless)  sky  (Wm"2);  CCLD  is  the  cloud  cover  fraction  (0-1.0);  Rsw  is  the  surface  reflection 
coefficient;  and  fSHD  is  the  shading  factor  of  vegetation  and  landscape  near  a  water  body.  If  the 
measurements  of  solar  radiation  are  available,  they  may  be  directly  applied  for  Jsw,  and  the 
calculations  of  Equation  (2)  may  be  avoided. 

Long  Wave  Atmospheric  Radiation,  Jlw 

J LW  =  ^  AIR®  ^ AIR  (1  +  0-  ^-7Ccld  )(1  ~  Rlw)~  £-wTR®TwTR  0) 

where  eAIR  and  ewm  are  the  emissivities  of  air  and  water,  respectively;  TMR  and  Twm  are  the  air 
and  water  surface  temperatures,  respectively;  R,w  is  the  surface  reflection  coefficient  of  long 
wave  radiation;  and  o  is  the  Stephan  Boltzman  constant  (5.670  x  10"8  W  m"2  K"4). 

Latent  Heat  Flux,  Je 


J F,  f  (UmD)(e AIR  eS  )  (4) 

where  f(UmD)  is  a  function  of  wind  speed,  UmD  ;  eAIR  is  the  air  vapor  pressure  (mb);  and  es  is 
the  saturation  vapor  pressure  (mb)  at  the  water  surface.  In  the  CMS,  f(UWND )  is  calculated  using 

/(^)  =  6.9 +  0.345^  (5) 
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where  UfmDJ  is  the  wind  speed  measured  at  7  meters  (m)  level  above  the  water  surface  (Edinger 
etal.  1974). 

Sensible  Heat  Flux,  Js 


(6) 


where  CB  is  the  Bowen  coefficient  (0.62  mb  K'1). 


Therefore,  the  total  heat  flux,  JT ,  in  Equation  (1)  can  be  obtained  by 


J T  J sw  J LW  ^  J E  "f”  J S 


(7) 


MODEL  ASSUMPTIONS:  CMS-Flow  is  presently  capable  of  2D  temperature  computations  in 
the  implicit  solver.  The  simulation  of  temperature  normally  requires  a  three-dimensional  solution 
due  to  the  presence  of  vertical  temperature  gradients  that  can  influence  the  flow  significantly.  It 
is  therefore  important  to  understand  the  limitations  of  2D  temperature  simulations  and  apply 
them  only  when  the  assumptions  inherent  in  2D  simulations  are  valid.  Typically,  2D  temperature 
simulations  are  valid  when  the  temperature  is  well  mixed  over  the  water  column.  These 
conditions  are  usually  met  for  shallow  bays  with  open  exchanges  to  the  ocean  or  gulf  and  strong 
tidal  signals  and  sufficient  wind  energy  to  provide  the  vertical  mixing.  Also,  the  assumption  of 
sufficient  energy  to  mix  over  the  water  column  is  valid  under  storm  conditions,  even  for  deeper 
water  bodies.  Finally,  when  the  exchange  with  the  open  sea  is  restricted  by  an  inlet,  the  tidal 
range  is  an  important  indicator  of  vertical  mixing  conditions.  For  low  tide  ranges,  significant 
vertical  stratification  can  occur,  even  in  shallow  bays  and  estuaries,  especially  when  the  winds 
are  calm.  Pritchard  (1955)  and  Cameron  and  Pritchard  (1963)  have  classified  estuaries  using 
stratification  and  freshwater  inflow  as  the  governing  criteria,  and  these  classifications  can  be 
used  for  guidance  in  applying  the  2D  simulations. 

TEMPERATURE  MODELING  IN  CORROTOMAN  RIVER,  VIRGINIA:  In  this  section, 
temperature  modeling  is  described  to  demonstrate  the  CMS  capability  in  the  Corrotoman  River, 
Virginia.  The  temperature  calculations  are  validated  against  measurements  at  two  monitoring 
stations. 

Background.  The  Corrotoman  River  is  a  small  tributary  of  the  Rappahannock  River  in  the 
coastal  plain  of  Virginia.  The  Y-shaped  river  has  two  branches,  the  western  and  the  eastern 
branch,  and  it  joins  the  Rappahannock  River  before  flowing  into  the  Chesapeake  Bay  (Figure  1). 
Freshwater  discharges  into  the  Corrotoman  River  come  from  a  few  small  creeks  and  the  annual 
flow  rate  is  very  small.  The  river  has  an  average  water  depth  of  1.7  m  with  a  mean  tidal  range  of 
0.35  m.  The  hydrodynamics  of  the  Corrotoman  River  is  predominated  by  tidal  process  with 
occasional  passages  of  meteorological  events.  Tide  and  wind  provide  sufficient  energy  to  mix 
water  vertically  in  the  shallow  river.  Observations  by  the  Virginia  Estuarine  and  Coastal 
Observing  System  indicate  that  depth-averaged  circulation  and  temperature  simulations  are 
applicable  to  the  system  as  the  temperature  is  well  mixed  over  the  water  column 
(http://web2.  vims,  edu/vecos/) . 
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Figure  1.  CMS  domain,  quadtree  grid,  and  bathymetry  of  Corrotoman 
River,  Virginia.  Triangles  indicate  the  survey  station  locations. 


In  application  of  the  CMS  to  the  Corrotoman  River,  a  quadtree  grid  system  was  developed  to 
discretize  the  computational  domain,  which  covers  approximately  13-16  kilometers  from  the 
river  mouth  to  the  upstream  of  the  western  and  the  eastern  branch.  The  downstream  boundary  of 
the  domain  reaches  to  the  7  m  depth.  Figure  1  shows  the  quadtree  grid  and  bathymetric  features. 
There  are  82,000  ocean  cells,  and  the  cell  size  ranges  from  20  m  in  the  upstream  to  80  m  in  the 
downstream  of  the  river,  which  permits  resolution  of  the  complicated  shoreline  and  narrow 
portions  of  the  river.  The  implicit  solver  of  the  CMS,  with  a  large  time-step  of  10  minutes,  was 
employed  for  simulations. 

The  temperature  data  used  in  the  study  were  downloaded  from  the  Virginia  Estuarine  and 
Coastal  Observing  System  (VECOS)  (http://web2. vims. edu/vecos A .  Two  sampling  stations  are 
located  in  the  Corrotoman  River  and  one  near  the  river  mouth  (Figure  1).  Continuous 
measurements  of  temperature  were  conducted  at  CRR004.02  with  an  interval  of  15  minutes,  and 
long-term  samples  were  measured  at  LE3.3  and  LE3.4  with  an  interval  of  3-5  weeks.  The 
temperature  data  inside  the  river  domain  (CRR004.02  and  LE3.3)  were  used  to  validate  the  CMS 
temperature  calculations;  data  near  the  river  mouth  (LE3.4)  provided  boundary  temperature 
conditions  along  with  the  WSE.  The  simulation  period  was  from  11  October  to  9  November 
2009.  Based  on  the  long-term  temperature  measurements  at  LE3.3  and  LE3.4,  a  constant  initial 
temperature  value  was  specified  for  the  entire  CMS  domain. 

CMS-Flow  was  driven  by  time-dependent  WSE  at  the  open  boundary,  wind  speed  and  direction, 
and  other  atmospheric  data  (solar  radiation,  air  temperature,  cloud  cover)  over  the  surface 
boundary.  WSE  forcing  was  downloaded  from  a  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  tidal  gage  at  Windmill  Point,  VA  (#8636580).  Wind,  air  temperature, 
and  cloud  cover  were  obtained  from  National  Centers  for  Environmental  Information 
(http://www. ncdc. noaa. sov/customer-support/online-store)  and  solar  radiation  data  from  National 
Solar  Radiation  Data  Base  (http://rredc. nrel sov/solar/old  dgta/nsrdb/)  at  Langley  AFB  station. 
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CMS  setup  for  temperature  calculation.  The  CMS  hydrodynamic  input  files  for  the 
Corrotoman  River  were  prepared  using  the  SMS  (Figure  2).  Instructions  for  setting  up  a  CMS 
temperature  simulation  are  described  below. 


Figure  2.  Files  for  the  CMS-Flow  temperature 
simulation. 


After  opening  “CR2.cmcards”  in  SMS  12.1,  right-click  the  simulation  name  under  “CMS-Flow 
Simulations”  in  the  SMS  Data  Tree  and  choose  Model  Control,  click  on  the  Salinity /Temperature 
tab,  and  check  the  box  next  to  Calculate  temperature.  A  default  time-step  equal  to  the 
hydrodynamic  time-step  will  be  specified.  In  this  case,  600  seconds  was  used  for  the  temperature 
calculation,  equal  to  the  10-minute  hydrodynamic  time-step.  Atmospheric  parameters  that  cover 
the  computational  domain  must  be  specified.  That  step  is  covered  in  the  Air-Water  interface 
section. 

1.  Temperature  initial  condition.  Because  of  the  large  spatial  variability  of  temperature  in  a 
coastal  system,  long  spin-up  periods  are  usually  required,  which  could  range  from  a  few 
days  to  weeks  for  a  simulation  to  reach  to  the  present  temperature  distribution.  To  shorten 
the  spin-up  time,  an  initial  temperature  field  should  be  specified  based  on  available  inform¬ 
ation.  There  are  two  options  to  assign  the  initial  temperature  condition  in  CMS-Flow: 

2.  A  global  initial  temperature.  Specify  a  constant  initial  value  for  the  entire  model  domain. 
The  temperature  value  can  be  specified  by  choosing  Constant  Water  Temperature  under 
the  Initial  condition  section  (Figure  3)  and  entering  a  value  under  Water  properties.  If  this 
option  is  applied,  it  is  best  to  define  an  average  representative  temperature  for  the  entire 
domain. 
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Temperature 

I?!  Calculate  Temperature 

Time  steps 

Transport  rate: 

600 

Initial  condition 


|  Constant  Water  Temperature  ▼  | 


seconds  ^ 


Water  properties 
Water  density: 

1025  kg/mA3 

Constant  Water  Temperature: 

13.35  deg  C 


Figure  3.  Temperature  dialog  showing  options 
for  the  Constant  Water  Temperature 
initial  condition. 

3.  Spatially  varying  initial  temperature.  Generate  a  spatially  varying  initial  temperature 
field  by  choosing  Spatially  varied  toggle  under  the  Initial  condition  section  (Figure  4a)  and 
then  the  Select  button.  Clicking  the  Create  Dataset  button  in  the  resulting  dialog  will  bring 
up  the  Data  Calculator  where  a  value  can  be  assigned  in  the  calculator  entry  box  and  then 
enter  an  output  dataset  name  such  as  “Initial  temperature,  C.”  The  Compute  button  must  be 
clicked  before  proceeding  out  of  the  dialog.  Close  all  the  way  out  of  the  dialogs  by  clicking 
Done  to  close  the  data  calculator;  then  click  OK  twice  to  close  the  CMS-FLOW  Model 
Control  window.  A  dataset  with  the  name  provided  in  the  calculator  will  appear  in  the 
CMS-Flow  data  tree,  as  shown  in  Figure  5.  Highlight  the  dataset  to  specify  different 
temperature  values  in  the  CMS  domain  in  the  same  way  to  modify  other  datasets  such  as 
“Manning’s  For  an  existing  dataset,  click  the  Select  button  under  the  Spatially  varied 
toggle  and  then  select  the  existing  dataset  for  the  initial  temperature  (Figures  4  and  5). 


Figure  4.  a)  Temperature  dialog  for  the  Spatially  varied  initial  condition,  b)  The  Data 

Calculator  for  specifying  the  value  and  the  dataset  name  of  the  initial  condition. 
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Figure  5.  SMS  data  tree  for  the 
newly  created  dataset. 


4.  Based  on  the  survey  data,  temperature  only  varies  by  0. 5-1.0  °C  from  near  the  mouth  to 
the  upstream  of  the  river  during  the  simulation  period.  Therefore,  the  initial  temperature 
was  obtained  by  the  temporal  interpolation  of  the  measurements  at  LE3.4,  and  a  constant 
value  of  18.85  °C  was  assigned  for  the  Corrotoman  River  system. 

5.  Temperature  boundary  conditions.  To  calculate  temperature  transport,  temperature 
values  at  CMS-Flow  boundaries  must  be  specified.  As  described  in  Li  et  al.  (2012),  there 
are  two  boundary  types  in  the  CMS  for  salinity  calculations:  WSE  boundaries  ( WSE - 
forcing  boundary)  (Figure  6)  and  freshwater  inflow  boundaries  (Flow  rate-forcing 
boundary).  For  temperature  calculations,  heat  transfer  at  the  air-water  interface  (surface 
boundary)  needs  to  be  calculated. 


Figure  6.  Temperature 

boundary  types  in 
the  CMS. 
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i)  WSE-forcing  boundary.  In  the  SMS  Data  Tree,  select  the  Boundary  Conditions 

coverage  in  the  Map  Module  and  use  the  Select  Feature  Arc  tool,  I,  to  left-click 
and  select  the  WSE  boundary  arc  (Figure  6).  After  selecting  the  arc,  right-click  and 
choose  the  Assign  Boundary  Conditions  option  (see  Figure  7a  for  the  Arc  Boundary 
Condition  dialog).  A  time  series  of  temperature  can  be  assigned  along  the  WSE- 
forcing  boundary  by  checking  the  Define  Temperature  Curve  box,  then  clicking  the 
Curve  undefined  button  at  the  bottom  center  of  the  dialog. 


Figure  7.  a)  Image  of  the  Arc  Boundary  Condition  dialog,  b)  Image  of  the 
temperature  boundary  input  window. 


6.  The  time  series  is  specified  either  by  clicking  the  Import  button  to  read  a  temperature 
boundary  input  file  in  xys  format  (Aquaveo  2010),  by  entering  time  and  temperature 
values  manually  in  two  separate  data  columns,  or  by  importing  temperature  data  from  an 
opened  Excel  file  using  Copy/Paste  (Figure  7b). 

7.  Figure  7  shows  the  temperature  measurements  at  the  long-term  station,  LE3.4  (September- 
December  2009),  which  were  assigned  to  the  offshore  boundary  with  the  WSE 
measurements  at  the  NOAA  tidal  gage,  Windmill  Point,  VA. 

i)  Flow  rate-forcing  boundary.  Following  the  same  steps  as  specifying  WSE-forcing 
boundary,  temperature  values  at  freshwater  inflow  boundaries  can  be  assigned  together 
with  flow  specifications. 

ii)  Air-water  interface.  Equations  (2)  through  (6)  are  used  to  calculate  heat  fluxes  at  the 
air-water  interface.  Air  temperature,  dew  point  temperature,  cloud  cover,  and  solar 
radiation  are  obtained  from  the  NOAA  weather  station,  Langley  Air  Force  Base, 
Virginia,  and  must  be  specified  in  the  Temperature  section  of  the  Salinity/Temperature 
Model  Control  interface  as  shown  in  Figure  8. 
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Atmospheric  Parameters 
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Figure  8.  Atmospheric  Parameters  dialog  in  Temperature  Model  Control  interface. 


SIMULATION  RESULTS:  For  the  demonstration  in  the  Corrotoman  River,  the  CMS-Flow 
simulation  was  conducted  for  a  30-day  period  (1 1  October-9  November  2009).  Depth-averaged 
current  and  temperature  fields  in  the  system  were  retrieved  from  four  snapshots  of  the  CMS 
results  at  a  6-hour  interval  on  21  October  2009,  corresponding  to  the  ebb  and  flood  currents, 
respectively  (Figure  9). 

Figure  9  shows  that  peak  currents  occur  near  the  river  mouth  with  a  current  speed  of 
1 5  centimeters  per  second  (cm/sec)  during  the  ebb  and  flood  periods.  Daily  variations  in 
temperature  are  clearly  displayed  in  the  figure,  which  are  more  significant  in  extreme  shallow 
(marsh)  areas  from  10  to  20  °C  and  less  than  1  °C  in  the  deep  portions  of  the  river. 

Model  temperature  calculations  in  the  system  were  checked  using  measurements  at  stations 
CRR004.02  and  LE3.3  in  Figure  10.  Stations  CRR004.02  and  LE3.3  have  water  depths  of  1.8 
and  6.6  m,  respectively.  Calculated  temperature  at  LE3.3  shows  smaller  temporal  variations  of 
less  than  5  °C  and  at  CRR004.02  a  larger  variability  of  6-7  °C  during  this  30-day  simulation 
period.  Controlled  by  meteorological  conditions,  the  temperature  at  the  shallow  site  also  shows 
large  diurnal  changes  of  more  than  2  °C. 

Statistical  parameters  were  calculated  for  the  evaluation  of  the  CMS  performance  at  CRR004.02. 
The  correlation  coefficients  and  the  root  mean  square  error  (RMSE)  are  0.866  and  0.523  °C, 
respectively,  which  indicate  a  good  agreement  between  the  calculations  and  the  measurements  at 
this  site.  The  sensitivity  tests  reveal  that  temperature  variations  in  the  system  are  primarily 
controlled  by  meteorological  conditions,  and  ocean  processes,  tides  in  this  case,  also  play  a  big 
role,  especially  in  deep  areas. 

The  CMS  simulations  represent  the  temperature  variations  in  the  Corrotoman  River  to  an 
accuracy  level  useful  for  understanding  general  temperature  patterns  in  the  river.  Based  on  the 
study  results,  further  improvement  can  be  achieved  through  increasing  knowledge  of  the 
temporal  variation  and  spatial  distribution  of  temperature  and  the  ocean  and  atmospheric 
processes  in  the  river/estuarine  system  (e.g.,  by  adding  a  term  to  account  for  the  flux  of  heat  into 
and  from  the  sediment  bed). 
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Figure  9.  Temperature  and  current  distributions  on  21  October 
2009  at  03:00,  09:00,  15:00,  and  21:00  EST. 


Figure  10.  Calculated  and  measured  temperature 
at  stations  CRR004.02  and  LE3.3  for 
October-November  2009. 
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CONCLUSIONS:  The  Coastal  Modeling  System’s  capability  in  computing  the  depth-averaged 
temperature  calculation  in  the  Corrotoman  River  was  demonstrated  in  this  technical  note.  The 
example  shows  that  the  procedures  to  set  up  a  CMS  temperature  calculation  are  straightforward 
and  user  friendly.  This  estuarine  application  provides  a  validation  example  for  the  CMS  and 
satisfactorily  reproduces  the  temporal  variations  and  spatial  distributions  of  temperature. 
Because  the  CMS  does  not  simulate  the  ice  freezing/melting  process,  the  temperature  calculation 
is  not  applicable  to  ice  covering  regions.  Further  improvement  to  the  CMS  can  be  made  as 
different  bulk  formulas  are  introduced  into  the  model  to  calculate  the  heat  flux  components  for 
applications  in  coastal,  estuarine,  or  riverine  environments. 

ADDITIONAL  INFORMATION:  This  CHETN  was  prepared  and  funded  under  the  Coastal  Inlet 
Research  Program  (CIRP)  and  was  written  by  Dr.  Honghai  Li  (Honghai.Li(a>nsace.  grmyjnil, 
voice:  601-634-2840;  fax:  601-634-3080)  of  the  U.S.  Army  Engineer  Research  and  Development 
Center  (ERDC),  Coastal  and  Hydraulics  Laboratory  (CHL).  The  CIRP  Program  Manager,  Dr.  Julie 
D.  Rosati  (Julie. D. RosatKmtsace. armv.mil),  and  Dr.  Weiming  Wu,  professor  at  Clarkson 
University  (wwuCcv, clarkson.  edu),  reviewed  this  CHETN.  Files  for  the  example  may  be  obtained  by 
contacting  the  author.  This  CHETN  should  be  cited  as  follows: 

Li,  H.,  and  M.  E.  Brown.  2017.  Temperature  calculations  in  the  Coastal  Modeling 
System.  ERDC/CHL  CHETN-IV-110.  Vicksburg,  MS:  U.S.  Army  Engineer 
Research  and  Development  Center,  h tty ://dx. doi. ore/ 1 0.21079/1 1681/21666 
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